Photonic structures with disorder immunity 
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Periodic and disordered media are known to possess different transport properties, either classi- 
cally or quantum-mechanically. This has been exhibited by effects such as Anderson localization in 
systems with disorder and the existence of photonic bandgaps in the periodic case. In this paper 
we analyze the transport properties of disordered waveguides with corners at very low frequencies, 
finding that the spectrum, conductance and wavefunctions are immune to disorder. Our waveguides 
are constructed by means of randomly oriented straight segments and connected by corners at right 
angles. Taking advantage of a trapping effect that manifests in the corner of a bent waveguide, we 
can show that a tight-binding approximation describes the system reasonably well for any degree 
of disorder. This provides a wide set of non-periodic geometries that preserve all the interesting 
transport properties of periodic media. 

PACS numbers: 42.25.Dd, 42.70.Qs, 73.63.Nm 



I. INTRODUCTION 

The study of transport properties in materials has im- 
portant appHcations in the technological realm. With the 
advent of metamatcrials fH, it has been demonstrated 
that such transport properties can be controlled by mod- 
ifying the structure of solids. Both compositional and 
geometrical parameters play a significant role in the de- 
sign of new materials. From the standpoint of physical 
phenomena, we distinguish the remarkable properties of 
photonic i2j|3| and phononic crystals jj], which emulate 
many features of electronic transport in solids: The pres- 
ence of bandgaps, the appearance of conical points in 
the frequency spectrum of multilayered structures and 
the realization of effective Dirac equations in lattices [5j 
are just a few examples. 

The aforementioned properties usually depend on 
the periodicity of the assembled structures, including 
quantum-mechanical realizations in one dimension. Such 
examples of tight-binding chains can be found in poly- 
mers and even in the modern construction of opti- 
cal lattices In this paper we go further and intro- 
duce disorder as another ingredient, with the purpose of 
constructing more flexible configurations. We show that 
the lowest energy band of waveguides with randomly ori- 
ented segments possesses the spectral properties of pe- 
riodic structures, such as bandgaps, conical points (or 
Dirac points) and a conductance band. These are clear 
indications of robustness under disorder. 

In connection with two-dimensional open systems 
and non-integrable billiards, we should mention that in 
some cases, randomly disposed boundaries and potentials 
give rise to wave-like manifestations of chaos |10|] in 
the form of level statistics [ll| among other signatures. 
Additionally, Anderson localization fl^ stands as one of 
the unmistakable signatures of disorder, affecting the cor- 
responding transport properties. 

In our study we establish a result in the opposite di- 
rection, namely that low energy waves in systems with 
corners are immune to the effects of disorder, with no 
localized modes in the lowest frequency band and a non- 



vanishing conductance band below threshold. In fact, 
the examples studied in this paper recover one of the 
paradigms of integrability and solvability: The homoge- 
neous tight-binding chain with nearest-neighbor interac- 
tions. 

The presence of corners in our disordered waveguides 
is of utmost importance. Chains of connected resonators 
can be proposed in many ways [islll^ . but the shape of 
such resonators and the number of supported resonances 
inside them has an important effect in the complete sys- 
tem, perhaps in a rather uncontrolled manner. The fact 
that a corner connecting two leads at a right angle al- 
lows only one bound state can be used to our favor. The 
trapping mechanism of a corner is of a purely wave-like 
nature and manifests itself at the lowest part of the spec- 
trum, where the wavelenghts are larger than the width of 
the guides. The presence of bound states in corners was 
first noted in [15] and their description was developed in 
[16<] by means of conformal maps and effective potentials. 

In this paper we study the spectrum, eigenfunctions 
and dimensionless conductance of disordered waveguides 
with corners. As an important result we obtain transmis- 
sion below the threshold of a straight waveguide (section 
III), forming a low frequency band located around the 
single trapped mode of an isolated corner or L-shaped 
waveguide and producing a gap which extends from the 
edge of the lowest band to the straight-guide propaga- 
tion threshold. Several realizations of two dimensional 
pipes are obtained by varying the orientations of their 
segments, showing thus the robustness and flexibility of 
the system. Then we give an analytical explanation of 
these effects by finding the effective interaction of a wave 
with a corner through a conformal map and proceed to 
connect corners in tigh-binding schemes (section IV). We 
finish with a summary and an outlook (section V). 



II. DEFINITION OF OUR SYSTEM 

We are interested in the transport properties of two- 
dimensional waveguides with corners bent at right angles 
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FIG. 1: (Color online) Wave amplitudes \'i>{x,y)\^ for influx 
coming from the left lead in waveguides (Model 1) with L = 
3d, N = 4, D = 5d. In (a) we set p = 0, while in (b) p = 0.5. 
In both cases the resonant energy E = 0.9412i?t was used. 



and randomly oriented segments. There are two models 
which can be considered. See Figs. [1] and [2] We describe 
their geometry as follows. 



1. Model 1 

This model is built by assembling blocks with the forms 
U and n. These blocks are randomly alternated as one 
moves along the array, with the only condition the ar- 
ray contains no loops. The resulting configurations are 
almost horizontal, minimizing the space. We relax this 
condition in Model 2. The parameters are the following: 
d is width of the waveguide, L is the length of the straight 
segments in units of d, N is the number of unit cells (n 
or U) forming the guide, the quantity p is the probability 
of finding U cells (for many realizations this is roughly 
the ratio of U to □ cells). 

Examples are shown in Fig. [TJ The waveguide is con- 
nected to two semi-infinite collinear leads of width D. 
The leads are attached to the waveguide by means of tri- 
angular contacts, avoiding strong difFractive effects. The 
width of the leads can be chosen arbitrarily, as long as 
waves of a wavelength larger than d can be supported. It 
suffices to take D = 5d. 



2. Model 2 

We consider waveguides with the same parameters as 
for Model 1, however here the cells forming the waveguide 
have forms L or F, allowing more flexible configurations. 
Now p is the probability of finding cells of type L. As 
before, we fix Z? = 5d. See waveguide examples in Fig. [2] 



A. Boundary value problem 




FIG. 2: (Color online) Wave amplitudes \'i/{x,y)\^ for influx 
coming from the left lead in waveguides (Model 2) with L — 
3d, N ^ 10, D = 5d. (a) p = 0, (b) p = 0.2, (c) p = 0.4, 
(d) p = 1. In all cases the resonant energy E — 0.9504_Et was 
used. 



tions which are defined by the shape of our waveguides. 
Our results shall equally hold in settings involving elec- 
tromagnetic waves [lo| or acoustic waves in the linear 
regime [l3l, as we are dealing with the Helmholtz equa- 
tion in a hollow guide. For a quantum particle of mass 
m we have 



2m " 



(1) 



and "0(2;, y) = at the boundary. In principle, we can 
infer the properties of our system by sending a wave of 
fixed energy through one of the leads (of width D) and 
finding the scattering amplitudes. The energy is given 
by 



E 



2m 



(2) 



where km and rrn: / D are, respectively, the longitudinal 
and transversal momentum components of the total wave 
vector with magnitude K — \/2mE / h. Our purpose is to 
explore the lowest energy region of the system. Therefore 
it is convenient to express the energy in the units E/Et, 
which normalize our quantities with respect to the energy 
of the lowest mode of the guide, i.e. the threshold energy 
Et = (?iV2m)(7rVd^). 



III. NUMERICAL RESULTS 

Using finite element methods we compute the scatter- 
ing matrix (^-matrix) which has the form 



We focus on the solutions of the stationary Schrodinger 
wave equation subjected to Dirichlet boundary condi- 



t r' 
r t' 



(3) 
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The symbols t, t' , r, and r' are M x M transmission and 
reflection matrices, where M is the highest mode given by 
the largest m beyond which the longitudinal wave vector 



becomes complex. Once the S'-matrix is known we cal- 
culate the dimensionless conductance (see p^| - |20| and its 
application to mesoscopic systems in Q)- 



T = Tr(tt^) 



(4) 



Results for Model 1 



In Fig. [3] we show the conductance T as a function of 
E for waveguides with N — 10, p — 0.5, and D — 5d. For 
comparison purposes we fix a single random sequence of 
unit cells and present results for L = 4.8c? and L = 5d. 
We note that the conductance plots change importantly 
for different values of L, but all the arrays display local 
maxima (resonant energies) in the region below thresh- 
old. The existence of a gap ranging from E ^ 0.95Et to 
E ^ Et is evident in Fig. ^a). 

In Fig. ^h) we have amplified the region of the con- 
ductance band. Here we observe an increase of the band- 
width as the geometric parameter L decreases, a behavior 
that can be interpreted as an increase of the coupling be- 
tween corners in a tight-binding regime. This shall be 
explained in further sections. The number of resonant 
peaks coincides with the number of corners in the array. 

In Fig. m we plot wave intensities |^'(x, y)^ at resonant 
energies for waveguides similar to the ones used in Fig. [S] 
There is a visible accumulation of wave intensities near 
the corners, becoming more pronounced as L increases. 
The accumulation, however, does not occur in every cor- 
ner: The intensities display the behavior of Bloch waves 
from corner to corner, despite the fact that the arrays 
are non-periodic, showing a quasi Id propagation in a 
periodic medium. 

We have found that the effect of disorder in the con- 
ductance band under scrutiny is minimal, as there are 
small differences in T for different configurations of the 
guides parameterized by the value of p. See section III.D. 



B. Results for Model 2 

As it is evdent form Fig. [21 the resulting configurations 
are more flexible than those of Model 1 (as long as the 
array does not intersect itself). This model allows to 
connect endpoints at arbitrary heights. The results for 
the conductance T as a function of E for waveguides with 
= 10 is quite similar to those previously described for 
Model 1. The resulting conduction band has resonant 
peaks that numerically approach T = 1 and formally 
reach T = 1. The number of peaks of T in the band for 
these conflgurations coincides with 2N. 



As before, we observe that once A'^ and L are fixed the 
conductance plots are almost the same for different val- 
ues of p. A more detailed description is given in section 
III.D. Now, the conductance band can also be analyzed 
as a function of L, increasing from 2.5d to 5d. A de- 
tailed analysis shows that (i) the bands are asymmetric 
for small L and the bandwidth A decreases with increas- 
ing L, see Fig.lSfa); and (ii) the band center E'ccntcr moves 
to smaller values of E with increasing L, see Fig. [S^b). 
See section III.C. 

From our calculations we can also conclude that the 
differences in the conductance profiles between ordered 
and disordered waveguides (i) are larger the smaller the 
value of L is; and (ii) are larger on the left side of the 
bands. 



C. Structure of the lowest conductance band 

Our goal is to study the low frequency bands formed 
by the conductance below threshold. To this end we an- 



H0.5 



(a) 1 


hi ' 1 




III — L = 4.8d 






1 mi 


III — L = 5d 1 


9 


0.95 1 



^ 10 




0.93 



0.94 



0.95 



E/E. 



FIG. 3: (Color online) (a) Conductance T as a function 
of E for waveguides (Model 1) with = 10, p = 0.5, and 
L — 4.8d, 5d. A band below threshold appears for both values 
of the intersite distance, (b) Amplification of (a) in the lowest 
conductance band. The bandwidth A increases as the corners 
approach each other, while the center -Bccntcr hardly moves 
when the geometry is altered. 



(d) 

FIG. 4: (Color online) Resonant wave amplitudes for influx coming from the left lead in waveguides (Model 1) with 

iV = 10, D = 5d, and p = 0.5. From top to bottom {L,E/Et) = (3d, 0.9412), (3. 2d, 0.9384), (3.4d, 0.9184), (3.6d, 0.944), and 
(3.8d, 0.9552). See Fig. E] 
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FIG. 5: (a) Bandwidth A and (b) band center -Bccntcr as 
functions of L/d, for waveguides of Model 2 with A'^ = 10, 
p = 0, and D = 5d. Both quantities decay exponentially with 
the distance. While the bandwidth A ranges from to 0.3i?t, 
the center varies only in the range 0.93i?t < -Eccntcr < 0.97Et. 



alyze numerically the variation of the bandwidth and the 
position of the center of the band as functions of the ge- 
ometry. Scale transformations of the guide in the form 
d — > Ad and L — > XL modify the spectrum trivially by 
an overall scale of A~'^. In order to modify the structure 
of the resulting spectrum we vary only one geometrical 
parameter, in our case L. In this way we find the be- 
havior of the band when the distance between corners 
increases. We take A as the bandwidth defined by the 
difference of the two energies below Et at which the quan- 
tity T < 10~^. The center of the band i^ccnter is simply 
taken as the energy which divides the number of peaks 
in two equal parts. See Fig. [3] 

In Fig. Owe find an exponential decrease of the band- 



width with the distance between corners. The val- 
ues range from (long distance interaction) to 0.3Et 
(strongest interaction at d = 2.5). This supports the idea 
that two coupled corners should be sufficnt to describe 
the coupling between sites and its behavior as a func- 
tion of the distance between them. On the other hand, 
the center of the band shows little variations: for long 
distances we have a center at 0.93i?t, while for d — 2.5 
the center approaches the threshold at 0.97i?(. It should 
be noted that at long distances, all levels tend to be de- 
generate at Q.93Et. The meaning of this energy shall be 
ellucidated in further sections, where we show that a sin- 
gle corner at a right angle can support only one bound 
state lying at such an energy. We shall use these features 
in order to establish the validity of a nearest-neighbor 
interaction in the theoretical description. 

D. Small effects at the edge of the band 

We have seen that the effects produced by disorder are 
small. Nevertheless such effects can be distinguished by 
a close inspection of the conductance at the edges of the 
bands. Here we analyze numerically the consequences of 
introducing disorder in the arrays. In Fig. we show 
the peaks forming the conductance band, with curves of 
different color for values of p = 0.1, ...,0.5. The effects 
due to disorder are not visible at this energy scale, even 
at a separation distance of L = 2.5d between corners. 
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FIG. 6: (Color online) (a) Conductance T as a function of E 
for waveguides (Model 2) with iV = 10, L = 2.5d, and D = 5d 
for p = (dashed curve) and 8 waveguide realizations with 
p — 0.1, 0.5 (color curves), (b) Deviation of resonant peaks 
of the conductance as a function of the number of defects. An 
increasing disorder produces a small deviation of the peaks at 
the edge of the conductance band. The largest shift is of the 
order 10~^-Et for p = 0.5 corresponding to the blue curve. 



In Fig. inUb), the scale has been increased, showing the 
effects on the position of the peaks at the lower edge 
of the conductance band. The shift of the peaks occurs 
downwards and it increases with the value of p, widening 
the conductance band by amounts less than 10~^£'f. 



IV. THEORETICAL DESCRIPTION 

Our aim is to describe the effects numerically found 
in the previous analysis. For this purpose, we follow 
three steps. 1) We describe the trapping mechanism of a 
single corner, ensuring the existence of one bound state 
below threshold. This shall be done by means of a con- 
formal map, giving rise to an effective interaction in the 
form of a position-dependent mass (in case of a quantum 
wire) or an effective dielectric function (for electromag- 
netic waves). 2) We analyze the interaction between two 
corners coupled in two different configurations, namely J" 
and n . As the spectrum of these two systems are approx- 



imately equal (even when distances between corners are 
reduced), we conclude that the differences can be treated 
perturbatively. 3) We introduce a nearest-neighbor tight- 
binding model where the atomic sites are represented by 
corners. The effect of disorder is described perturbativ- 
ley, explaining the effects at the edge of the conductance 
band. 



A. The trapping mechanism of a corner 

The L-shaped waveguide can be transformed into a 
straight one by means of a conformal map. The resulting 
Helmholtz equation acquires a position-dependent fac- 
tor in the Laplace operator. For quantum-mechanical 
waves this can be interpreted as a position-dependent 
mass, whereas for components of electromagnetic fields 
this can be thought as an effective dielectric medium. We 
proceed as follows. We solve the stationary Schrodinger 
(or Helmholtz) equation as a Dirichlet boundary value 
problem defined by our waveguide in Fig. [71 We use 
units 1 = h^/2m and coordinates a;, y to write 



I an 



0, 



(5) 



where is the interior of the array. In the following we 
describe our conformal map in order to find a position- 
dependent mass or an effective dielectric function of the 
coordinates. The obvious choice for a conformal set of 
coordinates is a function which maps an infinite straight 
strip into a bent waveguide. Let F be an analytic func- 
tion such that u = Ke[F{x + iy)],v — lm[F{x + iy)]. This 
leads to a transformation 



V2 



5(m,w)^2 



(6) 



d{x,y) 

where the Jacobian appears as a prefactor and satisfies 



d{u, v) 
d{x,y) 



dF{z) 



dz 



(7) 



(u.v) 



where z — x+iy. The boundary value problem ([5]) for the 
function ip{u,v) = (j){x{u,v),y{u^v)) has the equivalent 
forms 



il>{(),v)=ij{l,v) = 0, 
with an effective dielectric function r/(u, v) 
1 



(8) 



( .'iv+E i^(u,«) =0, 
^{0,v) = 1/^(1, v) = 0, 



(9) 
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with an effective mass fi{u,v) = . 

The passage from the wavefunction ip{u, v) to a nor- 
mahzable 4'{x, y) imphes the use of the Jacobian men- 
tioned above. Therefore, the old wavefunction 4'{x, y) 
would also satisfy the boundary conditions if the square 
root of the Jacobian does not contain strong singularities 
according to \imu^i^oip{u,v)\j^\(^u,v) = 0. The explicit 
form of F can be chosen in many ways. However its 
behavior near the corners is universal, since the angle 
formed by the walls of the array represents a branchcut 
of the map and it determines uniquely the rational power 
q appearing in F ^ {z — zq)'' , where zq is a vertex on the 
boundary. Here we find it convenient to construct F by 
means of a composition of two Schwarz-Christoffel trans- 
formations (see Fig. [T]). One of them maps an infinite 
strip of unit width to the semiplane and the other maps 
the semiplane to a tilted trigon [2l| . For waveguides bent 
in a right angle we have 



F{z) = u + iv = — arcsin 

TT 



37r\ 1 3 
z sm I — , — , — 

1 y ' 4' 4 



(10) 



where is the inverse of the regularized Beta function 
[2^ . For further details and more general bending angles, 
we refer the reader to [l6| . 

The spectrum of this system can be obtained by solving 
(HJ) or its equivalent forms ([U and It is composed by 
a single bound state lying below the propagation thresh- 
old of the straight segments and a continuum of energies 
above such a threshold. Henceforth we shall refer to the 
bound state energy and threshold energy as Eq and Ft 
respectively. It has been shown numerically and ana- 
lytically [2I, [m that Eo « 0.925^;* and that the wave 
function of the single bound state decays exponentially 
along the arms with a decay length A ^ a/So [13] • 



B. Two connected corners: J" vs □ pipes 

The results above suggest that each corner can act as 
an atomic orbital in a lattice, with the possibility of con- 
necting many corners to form a wide class of structures. 
In order to ensure the rise of a tight-binding chain, we 
focus first on the interaction between two corners in two 
different configurations: _r guides and □ guides. These 
two systems can be regarded as two-level atoms concern- 
ing their spectrum. They shall be studied numerically in 
order to show a level splitting (or coupling) which decays 
exponentially with the distance between corners. The 
degenerate levels approach the single bound state level 
of the L-guide, regardless of the orientation of the cou- 
plings (_r or n). The small differences in the properties 
of _r and n shall then be exploited in more complex as- 
semblies containing them as building blocks. In this way, 
the bandwidth and bandcenter of a disordered chain can 
be shown to be almost independent of the realization and 
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FIG. 7: (Color online) (a) Conformal coordinates 

obtained from (|10p in the form x(u,v) -\- iy{u,v) = 
csc(37r/8)/ (sin^ [(?i iu)7r/2] , 1/4, 3/4) . Its behavior near 
corners shows an abrupt change of the contour density, re- 
lated to the Jacobian of the transformation, (b) Contour plot 
of the Jacobian as a position-dependent mass, its value be- 
ing minimal around the position of the internal corner and 
approximately unity along the arms. 



can be further related to the interaction between two cor- 
ners (we can refer to them indistinctly as atomic orbitals 
or lattice sites). 

We solve the stationary Schrodinger (Helmholtz) equa- 
tion in two dimensions with the geometries J" and □ 
and Dirichlet boundary conditions. The finite element 
method finds all eigenfunctions and energies, particularly 
those lying at the lowest part of the spectrum. We show 
in Fig. |8] the resulting wave functions. The low-energy 
levels Fi and F2 , and the level splitting (or effective cou- 
pling) denoted by A = 2(i?2 ^ Fi), can be given as func- 
tions of the distance L between corners. The factor of 
2 has been introduced in order to compare A with our 
previous definition of bandwidth. The results are shown 
in Fig. O The exponential decay of the level splitting 
with the distance L ensures that the coupling of two L- 
waveguides emulates two-level atoms for both J" and □ 
guides. The small differences between the energies for 
_r and n configurations guarantee the immunity to dis- 
order of structures that are randomly built from these 
blocks. Given our numerical results, we can establish a 
model hamiltonian for energies exclusively below thresh- 
old. The two-level H has the form 



Hr = 



Eo 



A 

^ 4 



for the _r configuration and 



Hr 



Eo 



#(1 



#(1 + ^) 



£^0 



(11) 



(12) 
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(b) 




(c) 




(d) 




FIG. 8: (Color online) Bound states of two interacting cor- 
ners with L = 3d. (a) Symmetric state for X guide, (b) 
Antisymmetric state for J guide, (c) Symmetric state for n 
guide, (d) Antisymmetric state for n guide. 



for the n shape. The small differences in the level split- 
tings of the two configurations are given by Ae/2. The 
hamiltonian (fT2|) can be cast as a perturbation of pT|) 
i.e. Hn = Hj + eV with V an off-diagonal potential. We 
shall use this potential in the construction of chains with 
blocks of the _r and □ types. 
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FIG. 9: (Color online) (a) Level collapse as a function of 
the distance L between corners. The decay of the coupling 
A follows an exponential law. (b) Energy difference of the 
two low-energy levels for each configuration. Black curve: X 
guide, Red curve: □ guide. Discrepancies are almost negligi- 
ble, even near L/d = 0.5, where the difference between the 
splittings of the two systems is maximal i.e. Ae/2 ~ 0.017Et. 



C. The tight-binding chain with disorder 

Now that we have shown that the differences between 
the coupling of corners in X and n shapes are small, 
we can construct a disordered tight-binding chain with 
many of these blocks. In the following we shall use states 
localized at the corners as a basis for the Hilbert space of 
the wave operator for energies below threshold. Denoting 
the n-th corner state by \n) and its localized wavefunction 
by ^(2; — Xn) = {x\n), we propose the lowest energy band 
hamiltonian of a homogeneous chain of N corners as 

H^E„ + -J2{\n){n+l\+h.c.}, (13) 

n=l 

where only blocks of the X type (comprising two corners) 
appear. If N is sufficiently large, we obtain Bloch waves 
as eigenvectors of (IT3)) and the typical dispersion relation 
Ek = Eq + (A/2) cos fc for the corresponding eigenval- 
ues [131 ■ Now, in the presence of disorder introduced by 
blocks of the n type, we have the modification 



A 

=^0 + ^ ^{[l + efTpH]|«>(« + l|+h.c.},(14) 

where ap (n) takes the values and 1 randomly as a func- 
tion of the site n. The parameter p denotes the ratio of 
n to X blocks with < p < 1/2. We further impose the 
constraint that if crp(n) = 1 then (7p(n + 1) = in order 
to avoid self-intersection of the array. 

The resulting dispersion relation is depicted in Fig. (TU] 
for parameters Et = 1, A = 4 and e = 0.1. This set of 
parameters represents a strong perturbation in compari- 
son with the numerically obtained relation (Ae) / {2Et) — 
0.142 — 0.125 given by the difference between the black 
and the red curve in Fig. ^h). With this we show that 
the claimed robustness in a disordered waveguide is ac- 
tually stronger. Inspection of Fig. [TU] reveals that the 
bandwidth increases in small amounts with increasing 
disorder p. The states at the edge of the band suffer the 
greatest modification, whereas the presence of a conical 
point located at the center of the band is immune to dis- 
order. From the point of view of symmetry, the existence 
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FIG. 10: (Color online) Numerical energies as a function 
of the eigenvalue number n for a chain of 50 sites randomly 
coupled according to (|14p . The parameters are A = 4 and 
e = 0.1. The colors from red to blue indicate realizations 
corresponding to values of p increasing in steps of 0.1. The 
behavior of our disordered chain is consistent with numerical 
calculations of conductance peaks, as the bandwidth increases 
monotonically with the number of perturbations or defects, 
modifying the edge of the band. On the other hand, the 
conical point of the spectrum at eigenvalue 50 is immune to 
off-diagonal disorder due to the swapping symmetry. See the 
text below (fT7|) . 



of a conical point is protected by the fact that all cor- 
ners are equal, although they are not equally connected. 
Therefore, off-diagonal perturbations do not modify the 
swapping symmetry of corners. 



The energy shifts are given by 

A'-^^E,n = e{m\V\m) 

N-l 



Af-l 



-hO(l/7V3/2) 
= f^C0s(^)+0(W3/.)^ (16) 

where Np is the number of perturbed sites with n cou- 
pUngs, computed by summing over all ap{n). For a large 
number of reahzations, the ratio Np/N tends to p and 
we may write 



A^^'^Ejfi sa peAcos 



27rm 
N 



(17) 



Perturbative approach 



Interestingly, at the center of the band m — N/ 4 and the 
corrections vanish to leading order independently of p. 
Therefore, the conical points are protected. At any other 
region of the band, the shifts increase linearly with the 
disorder parameter p. The corrections become significant 
when (jl7p has a maximum, i. e. at the edge of the band 
corresponding to m = 0. The largest correction possible 
is therefore peA. This formula for the largest shift has 
an upper bound of 1.7 x lQ~^Et when p = 0.5 and with a 
maximum value of the perturbation Ae/2 = 0.017-E't at 
separation distance L = d/2. 

In order to discuss the predictive power of (jl7p . we 
can estimate the numerical energy shift of the blue peaks 
in Fig. [St^b) with respect to the dashed curve. For a 
separation parameter L — 2.5d, the difference of level 
splittings is of the order of 10~^ Et, as can be seen from 
the data plotted in Fig. IH^b). For maximal disorder we 
have again p = 0.5, leading to corrections at the edge of 
the band A^^^^Et = peA ^ 10~^ Et, in accordance with 
the numerical value. 



In order to describe the widening of the bands and the 
edge effects analytically, we find the corrections to the 
spectrum by using first order perturbation theory in e. 
To this end, we first approximate our unperturbed solu- 
tions by Bloch waves (finite size effects shall be treated as 
0{1/N) corrections). We have unperturbed eigenvectors 



m 



1 ^ 
1= exp 



r! = l 



yTTTD 

|„) + 0(l/iV), (15) 



where m is an integer. The Bloch quasi momentum is 
recovered in the limit A'^ — >■ oo in the form 27rm/A^ — >■ k. 



2. Wavef unctions and off-diagonal disorder 

When it comes to the discussion of localization, we 
should not ignore the fact that off-diagonal disorder 
in tight-binding arrays has been addressed before [25|, 
where a vanishing value of the transmission has been re- 
ported. The statistical transmission coefficient computed 
there, results in a decreasing function of the size of the 
system multiplied by the strength of some random poten- 
tial. Although our systems admit a wide range of config- 
urations, our random potential in (jl4p is perturbative. In 
the light of this result, localization lengths defined in con- 
nection with transmission coefficients are inversely pro- 
portional to our parameters e and p. As we have proven 
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that Ae/2 < 0.017£'f for all separation distances L, we 
can be sure that chains with thousands of corners do not 
feel the localization effects aforementioned. Moreover, we 
have used a dimensionless conductance in terms of the S 
matrix. Such quantity provides a more realistic approach 
to the transport properties, as it is directly connected to 
the numerical solutions of the scattering problem in a 
wire and it gives useful information for each realization 
of the system. 

V. SUMMARY AND OUTLOOK 

Let us summarize our results. We have studied numer- 
ically two models of waveguides with corners and seg- 
ments with random orientations. Immunity to such dis- 
order is a counterintuitive effect in the context of waves, 
and we verified its validity by showing that the lowest 
conductance band of the system is similar to that of a 
one-dimensional crystal. A theoretical description of the 
effects in question was given in terms of trapped modes 
in bent waveguides and their coupling in a tight-binding 
scheme. Finally, we showed that severely disordered con- 
figurations can be treated perturbatively at low energies. 
This enabled us to estimate the small effects at the edges 



of the conductance bands, while the immunity of coni- 
cal (or Dirac) points was confirmed at the center of such 
bands. 

Our study can be extended to other systems follow- 
ing similar principles. For example, upgrading to two- 
dimensional systems seems possible, as it only requires 
the existence of trapped states in cross-wires [Tl. We 
propose the emulation of atomic centers in a monolayer 
obtained by replicating the system periodically. In con- 
nection with particle statistics, we should mention that 
a canonical second quantization scheme can be proposed 
on our disordered one-dimensional lattice by promoting 
the localized or atomic states to field operators. Both 
fermionic and bosonic schemes are possible. This opens 
the possibility of describing the propagation of indepen- 
dent electrons in a quantum wire, charge carriers in a 
medium with an effective wave equation or photons of a 
fixed polarization and very low frequency. 
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